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Two general-relativistic hydrodynamical models are considered: a model 
of self-gravitating static configurations of perfect fluid and a model of steady 
accretion of fluid onto a black hole. We generalise analytic results obtained 
for the original polytropic versions of these models onto a wider class of 
barotropic equations of state. The knowledge about the polytropic solu- 
' tions is used to establish bounds on certain characteristic quantities ap- 

. pearing in both cases. 

"o 

1. Introduction 
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We will deal with two general-relativistic hydrodynamical models assum- 
ing the spherical symmetry: a model of self-gravitating static configurations 
of barotropic perfect fluids and a model of steady accretion of barotropic 
fluid onto a black hole. 

Although the first issue belongs to a classical astrophysics (see e.g. 
115) ). quite recently Karkowski and Malec gave analytic estimates onto the 
binding energy of the spherically symmetric, compact polytropes [7j. These 
estimates had been used to infer the connection between the maximal sound 
speed and the sign of the binding energy related to the stability of the 
configuration. 

The Newtonian version of the second model has been promoted to a 
general-relativistic one by Michel ,10; within the so-called test fluid approx- 
imation. The first fully general-relativistic description taking into account 
the backreaction (i.e., the impact the accreting fluid has on the space-time 
geometry) has been formulated by Malec [5]. 

Both Karkowski and Malec |7] and the authors of the works on the 
steady accretion cited above deal with the polytropic fluids, i.e., fluids with 
the equation of state expressed either as p = Kn T or in the form p = Kg r , 
where K and T are constant. Here p denotes the isotropic pressure, n 
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is a baryonic density, and g refers to the energy density appearing in the 
energy-momentum tensor of perfect fluid. 1 

In the present paper we show generalisations of these models obtained by 
considering barotropic equations of state, instead of the polytropes. More 
precisely: it follows that the knowledge of some properties of the poly- 
tropic solutions can be used to estimate those properties in a more general, 
barotropic case. 

Allowing for a wider class of equations of state constitutes a step toward 
more physical cases than just polytropes. A simple barotropic model of 
a white dwarf can serve as an example. The simplest white dwarf model 
is that described by a two component barotropic equation of state with 
the two components being essentially of the polytropic form: a polytrope 
with the exponent equal 5/3 for low values of the density and 4/3 in the 
ultra-relativistic regime. 

The order of this paper is as follows. In the next section we give ba- 
sic definitions of the quantities used in the presented analysis. Section |3] 
contains the discussion of the first model — that of static configurations of 
perfect fluids, where we deal with the estimates onto the binding energy in 
the barotropic case. Section |I] is devoted to the model of the steady accre- 
tion of barotropic fluids onto a black hole. At the end some closing remarks 
are made in section |21 



2. Basic equations and definitions 

In the present section we recollect some of the equations and definitions 
introduced in 0| and [7]. 

A general, spherically symmetric space-time can be described by a line 
element 

ds 2 = -N 2 dt 2 + adr 2 + R 2 (d9 2 + sin 2 d(p 2 ) , (1) 

where N, a and R are functions of a coordinate radius r and an asymptotic 
time variable t. 

We will consider a foliation of the space-time by the slices of a fixed 
time t. The only non- vanishing elements of the second fundamental form 
of such slices can be computed to yield K T r = dta/(2Na), Kg = Kt = 
dtR/(NR) (the convention for defining the second fundamental form used 
here is that of Wald ^B] ) • The trace of the second fundamental form reads 
trK = N^dt In {^JaR 2 ) . The analo gous calculations can be also performed 
for the two-spheres of constant radius r embedded in each of the time slices. 

We will refer to these two kinds of polytropes (coinciding in the Newtonian regime) 
as to a p - n and a p - g polytrope respectively. 



3 



The expression for the trace of the second fundamental form of such sphere 
is given by p = 2d r R/ (Ry/a). 

We assume further the energy-momentum tensor of perfect fluid 

= (p+g)u fM u u + pg» u , (2) 

where denotes the four velocity of the fluid, p is the pressure and g the 
energy density. 

The line element of the form (^Q) still has certain gauge freedom. While 
when dealing with static configurations one usually uses the so-called polar- 
gauge, i.e., one characterised by a condition trK = K£, for the description 
of the evolving fluids it proves to be convenient to introduce the comoving 
coordinates. This can be done by imposing an integral condition of the form 

1 \ 2 n /■OO 

-R(trK-K)j = — jf R' 2 {K?dR'+ 

1 f'OQ 1 /*00 

R' 2 KltiKdR' - — / R' 2 tiKdR'. 



2 -R Jr. 4i? j R 

One can check that this indeed gives u r = u e = = 0. 
The relation V^T^ U = leads to the continuity equation 

(guP) = -pV^. (3) 

Given a barotropic equation of state of the form p = p(g) it is, however, 
possible to find an integrability factor n (the so-called baryonic density) 
such that V« (nn^) = 0. It can be shown that for the energy-momentum 
tensor ((21) the function 

f dg' 

n = n\ exp / — — (4) 



g' + p(g') 



has the desired property. In turn, if we assume an equation of state ex- 
pressed as p = p(n) then the function g given by 

g = n + n[ n dn'^ ( 5 ) 
Jo n' z 

satisfies the continuity equation Q. Equations @ and © can be treated 
as expressing the relations between the baryonic and the energy density. 

For a given baryonic density n one defines the baryonic mass by a 
simple integration over the volume of the configuration. The baryonic 
mass enclosed in a ball B(r) of the coordinate radius r is thus equal to 
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^bar(^) = Jbm dVn = 8tt J q dR' (R'n/p). Here dV denotes the volume 
element induced on the hypersurfaces of constant time t. 

An analogous volume integral of the energy density defines the so-called 
rest mass M{r) = j B{r) dVg = 8vr j^ {r) dR' (R'g/p). Although it is clearly 
a geometric quantity on a fixed Cauchy hypersurface, it is not conserved 
during the evolution of the system, i.e., it can change between subsequent 
Cauchy time slices. 

A mass measure that, at least for compact systems, can be shown to be 
conserved (see e.g. jH]), and that is defined in terms of the energy density 

is the quasilocal mass m(r) = f B , r \ dV (Rpg/2) = Air J* dR'R' 2 g. 

The quasilocal mass of a static (and spherically symmetric) configu- 
ration of perfect fluid enclosed entirely within a sphere of an areal ra- 
dius R coincides with the asymptotic mass. It will be denoted here as 
ms = m(r(R)) = Air J Q dR'R' g. Similarly to iris, the symbols Ms and 
m bar w iH refer to the total rest mass and the total baryonic mass of the 
static, spherically symmetric system contained in a sphere of an areal ra- 
dius R. 

Let us also note that one can write down the line element for the static 
spacetime in the polar gauge as 

ds 2 = -N 2 dt 2 + — \^dr 2 + r 2 (d6 2 + sin 2 defy 2 ) . 

r 

In this case both N and m are functions of the (areal) radius r. 



3. Binding energy of the static perfect fluids 

The binding energy of the static, spherically symmetric and compact 
configuration of the perfect fluid is usually defined as a difference between 
its baryonic and asymptotic mass, i.e., -EWd = "Zbar — rns = JydVn — 
JydVg^/l - 2m/r. The paper by Karkowski and Malec [7] is mainly de- 
voted to give an analytic estimate of this energy for the polytropic configu- 
ration of fluid. However, one of the most useful results obtained there is a 
theorem bounding the potential energy — Uq defined as 



-U G = M s -m s 




dv- 2 " m 



1+ II- 2m 
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Theorem 1 (Karkowski, Malec) Consider a static spherically symmet- 
ric configuration of the self-gravitating fluid contained entirely within a 
sphere of the areal radius R. Let g ^ 0. Then 

-U G = M s -m s ^ 




Note that this is a very general statement — no explicit form of the 
equation of state is being assumed. 

For the sake of shortness we will denote by E the right-hand side of the 
above inequality. 

As an implication of this theorem one can write the following bound on 
the binding energy of the static, polytropic configuration 

E hind ^Ms-^-j-m s >Z-^j, (6) 

where P denotes the integral P = J v dVp. Below we will give the inequali- 
ties of this kind, valid for a wider class of barotropic equations of state. 

3.1. Binding energy of the barotropic fluid spheres 

In many cases considered in astrophysics barotropic equations of state 
do not differ much from the polytropes. This suggests writing a barotropic 
equation of state as 

p(n) =p(n ) ( — j , (7) 

where no denotes an arbitrarily chosen density of reference and the exponent 
r is a function of n. 

The function T introduced in this way appears to be a kind of a mean 
value of the so-called adiabatic index Fi = (n/p) dp/dn = ((g + p) jp) dp/dg. 
The two "gammas" are related by the formula 

n r^n') 



n f n 

r(n, no) = In -1 — / dn 

n Jna 



Clearly, every barotropic equation with n ^ and p ^ can be written 
in the form ((7J). We will, however, deal with a narrower class of equations 
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of state for which the function T is bounded, i.e., r m i n ?J T(n, no) ^ r max , 
and that the inequality r m j n > 1 holds. This implies that p(n = 0) = and 
also that p(n)/n — > as n —* 0. These properties will be used explicitly in 
further considerations. 

We present now two generalisations of © onto the case with the equa- 
tions of state of the form (J7J) . 

Lemma 1 Assume conditions of theorem Let the equation of state be 
of the form ^ and let the exponent T be nonincreasing with the baryonic 
density n. If the reference density no is such that n ^ no, then -Ebind ^ 



valid for n ^ no and T(n,no) being a nonincreasing function of n. The 
integral appearing on the left-hand side of the above formula can be easily 
computed to give 



and thus n ^ g — p(n)/(T(n,no) — 1) ^ g — p(n)j (T m i n — 1). It is now 
enough to use the bound on —U given in [7], i.e., fydVg — ms ^ H to 
obtain the desired estimate onto the binding energy, namely 



This completes the proof. 

The similar result can be formulated for the case with T(n, no) being a 
nondecreasing function of n. We decided to omit the proof. Although it is 
more complex than the preceding one, it can be based on the similar idea. 

Lemma 2 Let the assumptions of theorem Q hold and let T(n, no) be a 
nondecreasing function of n. Then 



a-p/ (r min - 1) . 



Proof. Let us start with the inequality 
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where denotes the Heaviside function and the A constant is given by 



Let us note that in lemma |2] no assumption is made on the value of the 
reference density no- 

Consider now the limit of no — > 0: (i) Since the volume of the region 
where n < no tends to zero so does the last integral occurring in (JSJ). (ii) In 
general by changing no we change the form of T(n, no). Suppose now that 
the expression r max — r mm remains bounded as no — ► 0. Remembering that 
we are still dealing with the equations of state for which p(n)/n — > for 
n — ► we see that the value of A appearing in (JHJ) tends to the unity, (iii) 
Assume further that T(n, no) is not only a nondecreasing function for a given 
no but is a nondecreasing function for all no — > so that the lemma |21 holds. 
Then in the above limit bound (JHJ) reduces to E'bind ^ S — P/ (r mm — 1) and 
the same inequality emerges as a result of both lemmas. 

The barotropic equation of state can be also expressed in terms of the 
energy density g. It is clear, that the barotropic equation of state expressed 
in the form p = p(n) can be transformed into the barotropic form in which 
the energy density is used, i.e., p = p(g); in some cases it is, however, useful 
to preserve the freedom of choice between these two dual descriptions. 

In the following we present a result, analogous to the preceding one, but 
formulated for the equation of state of the form 



As before we assume that the exponent T(g, go) is bounded, i.e., 1 < r m i n ^ 
T(q,Qo) s% T max for some constants r min and T max . 

Lemma 3 Let the conditions of theorem^hold. Let the equation of state be 
of the form and the exponent T(g, go) be a nonincreasing function of g. 
Let us also assume that the density of reference go is greater than the largest 
value of g within the whole configuration. Then #s E — P/ (r m j n — 1). 

The last statement presented in this section gives the suitable bound on 
the binding energy for the equations of state of type © with the exponents 
L(£, go) being nondecreasing functions of g. 

Lemma 4 Let the assumptions of theorem be satisfied. Let the equation 
of state be of the form (0) and let T(g, go) be nondecreasing function of g. 
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Under these conditions the following inequality holds 

A f 

Eund > AE + Am s -m s -- / dVQ(g - g )p + 

J- min — i JV 



r . 

A mm 

dV@(g - g)p{g ) 



v 



- 

\0o 



Here the constant A is given by 

j , | p(e>o)A _ ( r ™- 1 )( r m 1 ; n-i) 



£>0 

Similarly, under restrictions analogous to those that have been men- 
tioned in the comment to lemma 121 the go — > limit of the above inequality 
reads E hind > S - P/ (T min - 1). 

On the sound velocity and the sign of the binding energy 

The results presented in the previous section allow us to formulate condi- 
tions involving the sign of the binding energy that remain in a strict analogy 
to those given in [7] for the polytropic equations of state. 

Assume that the equation of state is given by ©• Then the sound 
velocity can be easily calculated to yield 

a 2 = pU(°-) dr ^eo) + r(g,g ) 



g J dg g 

, e \ dT(Q,Qo) p(eo) ( gY^- 1 

pin — + T{g,g ) 



go) dg go \go 

Let us introduce now the mean sound velocity a as 

P 



a 2 



Clearly, the quantity defined in this way is bounded from above and from 
below by the maximal and minimal sound speed respectively. Let us further 
define a? = PT m \ n /Ms- If the exponent T is a nonincreasing function of g 
and the density of reference go is chosen to be such that g ^ go or conversely, 
if dT/dg ^ and qq ^ g, then (In (g/go) dT(g, go)/dg + T(g, go)/o) 1 ^ 
Q/T{g,go)- Consequently j v dV (In (g/ g ) dT(g, g )/dg + T(g, g )/ g)' 1 ^ 
Ms/T m \ n . This, in turn, implies that a? a 2 . 
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Corollary 1 Under the assumptions described above, the inequality a 2 > 
r m in (r m in — 1) E/M5. is a necessary condition for the binding energy Sbind 
to be negative. 

Corollary 2 Assume that the conditions specified at the beginning of this 
section are satisfied. Then for the binding energy to be negative it is neces- 
sary that 

(CT \ r m i n — 1 
msj 

where a^ax denotes the maximum value of the sound speed. 

A brief discussion of inequalities of the type of inequality (|1U|) can be 
found in [7j. 



3.3. Physical examples 

So far we have only been dealing with fairly general assumptions about 
the equation of state. It is, however, easy to find explicit physical examples 
of equations of state that satisfy the assumptions we have made. Let us 
begin with nonincreasing "gammas". Consider a simple model of a white 
dwarf supported entirely by a degenerate electron pressure (see e.g. [TT]). 
The pressure of the noninteracting relativistic electron gas in the zero tem- 
perature limit can be described by the formula (we use the cgs units in this 
section to present the formulae in the most familiar form) 



4 5 
urate 

P- 



3h 3 



x{2x 2 - 3) VI + x 2 + 3 sinh" 1 
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where m e denotes the rest mass of the electron, c stands for the speed of 
light and h is the Planck constant. The parameter x appearing here is 
related to the number density of electrons n e by x = (m e c/h) (3n e /(87r)) 1//3 . 
For the white dwarfs that contain almost no hydrogen one can use the 
following relation between the number density of the electrons n e and the 
baryonic density: n e = n/ (2mn), where tbh denotes the rest mass of the 
hydrogen atom. In this case formula 1)11(1 defines a barotropic equation of 
state which can be trivially expressed in the form (J7| by taking T(n, no) = 
(lnp(n) — lnp(no)) / (Inn — In no). One can show that the function T is 
always decreasing; moreover, independently of the value of the density no 
we have T(n, no) — > 5/3 for n — * and T(n, uq) — * 4/3 for n — * oo. 

Examples of the equation of state with nondecreasing barotropic indices, 
satisfying conditions of statement [21 can be found in ^1 or [5] . These 
equations had been used for the description of the stiffening of matter at 
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supranuclear densities. They are essentially of the form (JJJ) with the ex- 
ponent T given by T = r m i n + r]Q(n — n nuc ) In (n/n nuc ), where r m ; n and rj 
are some constants, n nuc denotes the nuclear density and G is the Heaviside 
function. The density of reference no is chosen to be equal unity in the used 
system of units. Let us note that the function T defined in this way will 
still be nondecreasing, even if we allow the reference density no to tend to 
zero. Additionally, if we assume r max = r(n max ,no) for some fixed n max , 
then r max — > r m j n as no — ► 0. Thus, this kind of equation not only fulfills 
the assumptions of lemma El but it also assures the desired limit behaviour 
at no — ► 0. 

4. Quasistationary (steady) accretion 

We will consider a spherically symmetric cloud of gas falling onto a non- 
rotating black hole. The black hole provides the simplest choice for the 
central object as one can assume that no shock waves occur. 

In this chapter we will work in comoving coordinates. Let us define a 
fluid velocity U by U = d t R/N = R (trK - K r r ) /2. This allows us to write 
the Einstein equations in the form (see e.g. [7j) 

pR=2 fZ^m^ (12) 

(the Hamiltonian constraint) and 

d R (R 2 U) - R 2 tiK = (13) 
(the momentum constraint). The evolution equation reads 

d t U = \( P R?d R N - - AnRNp. (14) 

The continuity equation © can be now written as 

d tQ = -NtiK(g + p). (15) 

The conservation of the energy-momentum tensor provides us also with the 
relativistic version of the Euler equation 

Nd R p+(g + p)d R N = 0. (16) 

We will assume that the accretion is steady and relatively slow, so that 
the central mass can be regarded as constant (see e.g. |3]). More precisely: 
(i) the accretion rate, defined as rh = (pt — (dtR)d R )m(R) for the given areal 
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radius R, is assumed to be constant in time; (ii) the fluid velocity U, energy 
density g, sound velocity a etc. should remain constant on the surface of 
fixed R: (dt — (dtR)d R )X = 0, where X = U, g,a, . . . 

The accretion rate introduced above can be computed to yield 

m(R) = -AirNR 2 U(g + p). (17) 

It was proved by Malec [S] that under the above conditions the accretion 
rate is independent of the surface (characterised by a given R) for which it 
is calculated, i.e., d R m = (cf. equation (|23]> below). 

Under the above assumptions the original partial differential equations 
become ordinary ones. The evolution equation (|14|) can be written as 
NUd R U = (pR) 2 d R N/A - Nm(R)/R 2 - AnNp, while differentiation of 
equation JEJ) gives NUd R U = N(Rp)d R (Rp)/4 - Nm(R)/R 2 + inRNg. 
By subtracting these equations we obtain one of the key equations of the 
presented model, namely 

Finally, under the assumption of the quasistationarity of the accretion, the 
momentum constraint ()1H|) together with the continuity equation yields 

d R ln(\U\R 2 ) = -^. (19) 
' g + p 

Both of these equations can be integrated starting from the outer boundary 
of the accreting cloud. If iVoo, p^ and Roo denote the adequate values at 
the boundary, we can write 



pR PooRoo V Jr v R' 

Similarly integrating equation ()19|) we get 



U =l^i ( 21 ) 
R z n 

where A is an integration constant. Analogous algebraic form can be given 
also to the Euler equation ((TH|) . namely 

N = ^. (22) 
g + p 

Here B denotes another integration constant. One can check the validity of 
this formula using Q). 
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Let us note that formula (|T7|) reduces to 

m = -4vrAB. (23) 

In principle, one can use rh instead of a constant A. This is what was done 
in the classical analysis by Bondi [Ij; we will, however, leave the constant 
A as it is, in order to preserve the clarity of the equations. 

In what follows we will also put = 1. This can always be done 
provided that a suitable foliation of the external Schwarzschild space-time 
is chosen. 

Let us now define a sonic point as such, where the length of the spatial 
velocity vector equals the speed of sound \U\ = 2U/(pR) = a. In the New- 
tonian limit the above definition coincides with the standard requirement of 
the equality between the velocity of the fluid and the local sound speed. In 
the following we will denote by the asterisk all values referring to the sonic 
point. 

Below we show a theorem being a simple generalisation (onto the barotropic 
fluids case) of Malec's result for the type ft - Q polytrope. 

Theorem 2 For the barotropic fluid accreting in the quasistationary way, 
the following equation holds at the sonic point 

°-( 1 "^ +<! -) = c, * 2 = ! w +c " (24) 

where c* abbreviates the expression c* = 2irR 2 p if . 

Proof. The equation Q18|) can be easily written in the form 

d R \nN = (~pR 2 d R (Rp) + 4itR 2 q + AirR 2 ft 

Now, differentiation of the squared Hamiltonian constraint equation (fT2 
gives 

]pR 2 d R (Rp) +4ttqR 2 = ^&- + RUd R U = 
4 R 



and thus 



On In N = -±3 - 2U 2 + ^ d R (U 2 R*) + 4^p) 
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An other expression for the logarithm of the lapse derivative can obtained 
from the Euler equation ()16|) . namely BrIiiN = —dnp(g + p). Introducing 
the sound velocity and remembering equation (|19|) we get (under assumption 
that the fluid is ruled by the barotropic equation of state) Or (U 2 R 4 ) = 
(2/a 2 ) 8rN. The above equations lead easily to a key relation 

(, ^U 2 \_ /rT2 „4N IGU 2 R fm(R) o n -\ 
('-iwj 88 (U R ) _ - jp- ^ - V> + 2«R^ . 

This equation was first obtained by Malec |Hj , but with an explicit use of the 
polytropic equation of state. The proof can now be completed immediately 
by writing the above equation together with the equation (j!2j) at the sonic 
point. 

4-1. A case without backreaction 

If the mass of accreting fluid is negligible compared to the total asymp- 
totic mass m, i.e., 

4vr / dR'R' 2 S <^m, (25) 

J R>2m 

equation (|2U|) yields N rj pR/2 « y*! — 2m / R + U 2 . Now, under an addi- 
tional condition c* = 2ttR 2 p^ <C 2m/ R* equation (|24jl can be rewritten in 
the form al (1 - 3m/ (2i?J) = U 2 = m/ (2JS*). Notice, that due to <dSJ) we 
have m(R) ~ m for i? > 2m and we are, in fact, considering the motion of 
the infalling fluid in the fixed Schwarzschild space-time (the so-called test 
fluid approximation). 

The steady accretion model of Malec has its continuation in a recent 
paper by Karkowski et al. [H] . The main aim of this work was to investigate 
the effects of backreaction of fluid onto the metric. It appeared that, at least 
in the polytropic case (p - g polytrope), some parameters of the sonic point 
(a 2 , U* and the ratio m^/R*) are the same in the full model with back- 
reaction as in the test fluid approximation, provided that the asymptotic 
parameters are the same in both cases. This shows that for some purposes 
the knowledge of much simpler test fluid approximation is satisfactory. Al- 
though this fact has been proved only for the p - g polytropic equation of 
state, the numerical calculations show that it holds true also for the p - 
n polytropes. We conjecture that it is also valid for a wider subclass of 
barotropic models. 

We will deal with the test fluid approximation until the end of this 
paper. Among the presented results the special emphasis will be given to 
the bounds on the speed of sound at the sonic point for a class of barotropic 
equations of state. If the above conjecture is true, this result will hold also 
for the general case that includes the backreaction. 
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4.1.1. The existence of the solutions 

Within the test fluid approximation the following existence theorem can 
be proved. 

Theorem 3 Assume that there exists a sonic point, the backreaction effects 
can be neglected, the equation of state is a barotrope given by Q) with the 
restriction that 1 < T(n, no). If in addition the inequality 

n-^-lna 2 < 2 f a 2 + - | (26) 
dn \ 3 / 

holds at the sonic point, then outside the event horizon there exist at least 
two solutions that bifurcate from the sonic point. 

Proof. Let us start with the Euler equation 1)22(1 expressed in the form 
Ls i±l = l = B sP . ( 27) 

n N /, 2m i A 2 

Y r~ 

The functions L and P denoting the left and right-hand side of this equation 
play a crucial role in the entire proof. 
From relation © we have 

L = l + 1 + [ n dn>^. (28) 
n Jo n' 2 

If the equation of state will be of the form ((7| and the restriction 1 < 
r m in ^ r(re, no) will hold, then, as it has been already stated, p/n —* as 
n — > 0. Additionally, the integral appearing in the formula (|28|) converges 
and lim n ^o J Q n dn'p(n')/n' = 0. Due to these facts we conclude that L —> 1 
as n — > and, what follows simply from the equation (J27|). P — > 0. The other 
limit, i.e., that of n — > oo gives, basing on the same assumptions, L — > cxd 
while for P we have lim n ^oo P = B/^Jl — 2m/R < oo. Thus in both cases, 
n — > and n — > cxd, the same inequality L > P holds independently of the 
value of R > 2m. Now, if it happens that for some value of n (denoted 
as n x ) and a fixed value of R an opposite inequality L < P is satisfied, 
it follows from the continuity argument that the graphs of L and P must 
intersect at least in two points, i.e., at least two solutions exist at R. 

The suitable value of n x for a given R can be chosen as n x = n* {R/R*) x . 
For the function n x (R) one can show that at R = R* we have dnL{n x ) = 
dnP{n x ). Careful computation of the second derivatives gives in turn 



2 



d 2 R P(n x )\ R=Rt = P(R*)j^ {{Sal -2)x 2 -9x-6 
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and 

,2 



<^M R=i?t = L(R^ ( ( al + (n ^Ina 2 ) 



R=R, 



It is easy to observe that the sign of the expression d R (P(n x ) — L(n x )) taken 
at the sonic point is equal to the sign of 

2al - 2 - ( ni- In a 2 ] | x 2 - S.r - G. 



dn / R=R* 

The discriminant of the above polynomial of x is always positive provided 
that inequality (|2G|) holds at the sonic point. In this case it is possible to 
choose such value of x that d R L(n x )\ R=R < d R P(n x )\ R=R . This in turn 
ensures that at least in some vicinity of i?* we have L(n x ) < P(n x ). 

The domain of the existence of the solution can be extended to the whole 
domain outside the horizon by exploiting the observation that 

9 „ (i _ P) = ^( 1 _^). (29) 

The argument allowing us to conclude about this extension is as follows. 
Assume that there exists R\ (say R\ > R) such that both solutions exist 
for R* < R < R\ and that there is no intersection of the graphs of L[n) 
and P(n) for R > R\. Again it follows simply from the continuity that for 
R = R\ there exist just one point where the two graphs intersect. Obviously, 
the derivatives with respect to n of both L and P would be equal in this 
point. That, in turn, contradicts formula ()29j) giving a non-zero value for 
d n (L — P) as long as R / i?* and as long as we are dealing with a finite 
baryonic density n and a positive value of the sound speed. Since we have 
already shown that for n — > oo we have L ^ P in the region outside the 
horizon, the baryonic density must remain finite and we are left with the 
conclusion that at least two solutions exist in the domain described by 
R£ (2m, R^). 

It should be stated that for both polytropic equations of state, i.e., for a 
type p - g and for a type p - n polytrope condition (|2Gjl is satisfied provided 
that the suitable polytropic exponent is less than 5/3. 



4.1.2. The polytropic equation of state 

The original model of Malec |H] uses the polytropic equation of state 
of the p ~ g type. We will adhere now the alternative and probably more 
popular equation of state 



p= Kn r , 1 < T ^ 5/3. 



(30) 
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Introducing a new function £ denned as £ = TKn r ~ 1 /(T — 1) one can 
immediately integrate the Euler equation (|T5|l to obtain N(l + £) = 1 + S,oo, 
where £oo denotes the value of £ at the outer boundary of the cloud. 

The sound velocity can be also conveniently expressed with the use of 
£, a 2 = (r - 1)£/(1 + £). It follows that 

a 2 < r - 1 (31) 

is a strict bound for the sound velocity No such bound exists for the 
alternative p - g polytrope. 

The Euler equation integrated above can be also expressed in terms of 
the sound velocity 

r _ i _ a 2 = N ( r _ i _ a 2 o ) . (32) 

Notice that the square of the speed of sound is directly proportional to the 
lapse function N here. This fact results with the significant simplification 
of the equations when compared to the p - g model. 

It has been already shown in ;2] that there always exists a unique sonic 
point for the p - g model in which the backreaction is neglected. One can 
parallel Malec's proof also for the p — n case simply by analysing the sonic 
point equation 

r - 1 - al = iV* (r - 1 - c&) . (33) 

It follows that there always exists a unique sonic point in the p - n model 
without the backreaction (i.e., a unique solution to the above equation) 
provided that V 5/3. Moreover, it must be located outside the horizon, 
i.e., R* > 2m. 

In this point the p - n model differs from the p - g model. While in 
the p - g model the sonic point may exist, at least formally, below the 
horizon 2 , the p - n equation of state does not allow (even theoretically) 
such a possibility to happen. 

Further progress can be made by considering equation (|33|) squared and 
rewritten in the form 

3 {alf + (7 - 6r) (al) 2 + (T - 1)(3T - 5)a 2 + 
+«L (2(r - 1) - al) = 0. (34) 

2 Such situation could contradict the established views on the properties of matter. If 
the square of the sound speed in the p - q polytropic model exceeds the polytropic 
index Y we have p > g, i.e., one of the energetic conditions breaks down. This is since 
for p - Q polytrope one has 
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This is a purely algebraic cubic equation for a*. It is easy to show that 
for T ^ 5/3 and for the interesting range of a 2 ^, i.e., < a 2 ^ < T — 1 three 
real different roots of this equation exist. Since the original unsquared 
equation ()33j) can be proved to have a unique solution, only one of the 
mentioned roots of (|34j) is the one we search for. Careful analysis shows 
that the first of two remaining solutions lays outside the range < al < 1 
and the second simply does not satisfy the original equation (|33|) . The only 
physical solution can be expressed by the following Cardano's formula (see 
also @]) 



al = i{ (>r - 7 + 2(31' - 2) cos 



7T 1 f 1 / 

- + - arccos^ — — ^ 54T 6 + 

3 3 \2(3r-2) 31 



-351P 4 - 558r + 486(r - l)a£, - 243a^ - 259 



(35) 



One may check, that for a 2 ^ = T— 1 the above gives al = a 2 ^ and for a 2 ^ = 
we get al = 0. 

The formula (|21[) for the value of the velocity U at a given radius R 
can be transformed into the form involving the sound speed instead of the 
baryonic density. For the p — n polytrope we have 



U = U* 



Rl falT-l-a 2 \ ^ 



R 2 \a?T — 1 — al 



This equation, together with the Euler equation (|32j) and equations (|2Uj) 
and (|12|). constitutes the complete set describing the full polytropic model 
that takes into account the backreaction effects. Since in the test fluid 
approximation we have obtained an analytic expression for the sound speed 
in the sonic point and 



a 



I „r 1 ^ ml + 3al 



U *=T^^^ N *= r^-^ , R* = ^^i^, (36) 

at 



l + 3a2' ^1 + 3a2 ' 



i.e., we know precisely all quantities characterising the sonic point, it is pos- 
sible to express the accretion rate with just one analytic formula depending 
on asymptotic parameters a^ and Uoq. Indeed, after some calculations one 
can arrive at the form 



m = -4irNR 2 U(g + p) = 
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In this point the model with the p - n polytrope appears to be simpler than 
the alternative p ~ g model in which the value of the accretion rate can be 
only estimated (or computed numerically) even when the backreaction is 
neglected. The detailed comparison of the two models can be found in [S|. 

Let us also notice that some authors (see e.g. ^D) define an accretion 
rate as the time derivative of the baryonic mass (instead of the quasilocal 
one). The accretion rate defined this way can be expressed by the for- 
mula rhbar = (dt — (dtR)dR)mY )ai (R) = —8irNURn/p. In the test fluid 
approximation where N ~ pR/2 we get TObar = —^R 2 Un. Finally, equa- 
tion 1)22(1 allows one to show the following relation between the two rates 
rh = Bulbar = ((goo + Poo )/ n oo) w-bar- Both these accretion rates become 
equal for g^ = n^, p^ <C g^o, i.e., when the fluid at the outer boundary is 
nonrelativistic. 



4.1.3. Estimates onto the speed of sound 

The analytic solution of the sonic point equation for the p — n polytropic 
model in the test fluid approximation allows one to make estimates onto the 
sound speed at the sonic point for a wider class of barotropes, namely for 
those with monotonic barotropic exponents. 

Indeed, let us assume a barotropic equation of state expressed as 

, r(n,n_) / \ r(n,n+) 



p(n)=p(n-)\^—\ =p( n+ )^—j , (38) 

where densities of reference n_ and n+ are chosen in a way which ensures 

that n_ ^ n and n + n (i.e., n_ should be less or equal to the smallest 

baryonic density appearing in the model 3 while n + should be greater than 

the largest baryonic density 4 ). Let us additionally assume that 1 < r min ^ 

r(n,n_) ^ r" ax and 1 < T+ in < T(n,n+) < T+ ax for some constants r~ in , 
r>— r+ aTlr | p+ 

x max) A min allKX ^ max' 



Lemma 5 Consider a spherically symmetric, steady accretion onto the black 
hole with a negligible backreaction. Let the equation of state be expressed in 
the form h°3$\) and let both exponents T(n, n_) and T(n, n + ) be nonincreasing 
functions of the baryonic density n. Then the square of the sound speed at 
the sonic point is bounded by two analytic expressions a 2 _ ^ a 2 ^ a\, where 



.2 1 



a_ 



6r min -7 + 2(3r+ n -2) cos 



/ , 

3 + (39) 



3 This means, in practice, n_ ^ rioo 

4 For most cases it should suffice to have n+ n(2m) 
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arccos < 



2 (3r: 



2)' 



2(27(r+ n ) 3 -54(r+ n ) 2 + 



+■ 36r+ n - 8) - 243 (r min - i) 2 ( r ™_ 1 a °° 

i max 1 



2 \ 2 



and 



i or; 



+- arccos < 



+ 36r; 



7 + 2 (3r max 
1 

2 (3r max 2) 

8) - 243 (T- 



2) 



cos 



7T 

3 + 



(40) 



2 27fr: 



1 



-54 (r max ) 2 + 

l-a; 



1 



A completely analogous lemma can be formulated for nondecreasing 
barotropic exponents. 

Lemma 6 Let the conditions of lemma\^be satisfied with an exception that 
this time both exponents T(n, n_) and T(n,n+) are assumed to be nonde- 
creasing functions of n. Then the square of the sound speed computed at the 
sonic point is bounded from above and from below by analogous expressions 
a 2 ^ a 2 ^ a\, where 



. 6r mi„-7 + 2(3r min -2) cos 



+ 



1 f 1 If 

- arccos < 5- 2 

3 \2(3r- in -2) 3 V V 

- 8) - 243 (r min - 1 



7T 

3 + 



3 

+ 36r 



27fr min ) 3 -54(r min ) 2 + 



r+ -l-a 2 ^ 2 ' 



and 



I or; 



+- arccos < 



7 + 2(3r max -2 cos 



7T 

3 + 



2 3T1 



2)' 



2(27(r max ) 3 -54 (r max ) 2 + 



+ (36r max - 8) - 243 (r max -i)' 



l-a 



^min ^ 
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The proofs of these lemmas have been omitted since they are straight- 
forward. 

4.1.4. Bounds on the accretion rate in the barotropic case 

Having estimated the sound speed at the sonic point we should be able 
to obtain bounds on the accretion rate m as well. We will restrict ourselves 
to the case of lemmaOD The bounds on the accretion rate for the alternative 
case described in statement El can be obtained in an analogous way. 

Equations (|36|) (valid in the test fluid case) give the values of £/*, iV* 
and R* expressed in terms of the sound speed a*. Thus, as it is clear from 
formula (|17j). the only quantity that needs to be estimated is 

Suppose now that the conditions of statement \5\ hold. One can find 
T(n,n + )p/a 2 ^ g + p ^ Y(n,n-)p/a 2 . Furthermore, making use of the 
same inequalities together with their asymptotic versions (i.e., at R = Roo) 
we can write 



where no stands for re+ or 

It follows that the bounds on the accretion rate rh can be formulated as 




r ("."o) 



and similarly 





r 



!;:/!. \ 



< irm 



2 l + 3a^ T 
a* < 



iax ~i \ I n oo 

—P{n-) 

4 \ n - 



) 




a: 



Vl + 3ag r 




max 



) 



r. 



max 



-1 



and 
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Obviously, the above formulae together with formulae (|4()j) and (|41jl giving 
the lower and upper bounds for the value of al complete our task. One 
can also check that these inequalities are exact. They are saturated for the 
p - n polytropes in which case their right-hand sides reduce to the same 
expression (|37|) . 

5. Closing remarks 

The stability of the self-gravitating polytropes against radial perturba- 
tions has been examined by Tooper |15j . All configurations with the neg- 
ative binding energy happen to be unstable, however the positivity of the 
binding energy is only a necessary condition of the stability of the configu- 
ration. The problem of the stability of the self-gravitating configurations of 
general barotropic fluids remains (to our knowledge) still open; it is, how- 
ever, intuitively tempting to treat the sign of the binding energy of such 
configuration as carrying information about the stability. Thus inequalities 
presented here, relating the maximal speed of sound within the considered 
configuration and the sign of its binding energy, might be treated as some 
versions of the Jeans inequality. These inequalities, being straightforward 
generalisations of those given by Karkowski and Malec constitute the nec- 
essary conditions for the negativity of the binding energy. Generalisation of 
the sufficient condition formulated in [7j remains still an open task. 

As to the steady spherically symmetric accretion we believe that at 
least from the relativist's point of view the issue of outermost interest is 
the backreaction of the fluid onto the space-time metric. This effect has 
been investigated so far only for the polytropic case UJ. In this paper we 
were mostly dealing with the test fluid approximation obtaining information 
about the characteristics of the sonic point for a class of barotropic equations 
of state. Since the characteristics of the sonic point are expected to be the 
same in the full picture that includes the backreaction, our results should 
be relevant in the general, spherically symmetric case. 

One of the most striking effects caused by the backreaction is that the 
accretion rate achieves a maximum value for the ratio of the mass of the 
fluid versus the total mass of the system being equal to 1/3. This result 
has been obtained both for the p - g and p - n polytrope. We conjecture it 
holds true also for a wider class of barotropic equations of state. 
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